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, Abstract 

The Turaev-Viro invariants are scalar topological invariants of three-dimensional manifolds. Here we 

■ show that the problem of estimating the Fibonacci version of the Turaev-Viro invariant of a mapping 

^5 ' torus is a complete problem for the one clean qubit complexity class (DQCl). This complements a 

, previous result showing that estimating the Turaev-Viro invariant for arbitrary manifolds presented as 

T-H ■ Heegaard splittings is a complete problem for the standard quantum computation model (BQP). We also 

CO ' discuss a beautiful analogy between these results and previously known results on the computational 

, complexity of approximating the Jones Polynomial. 

S: 

^! 1 Introduction 
> , 

c : . . 

, Classifying the power of quantum computers is a fundamental problem in quantum information science. The 

^ ■ computational power of a general-purpose quantum computer is identified with the complexity class BQP 

^3^1 (bounded-error quantum polynomial time). The famous problems of factoring and discrete logarithm, for 

instance, are in BQP. An essential ingredient of BQP computation is the ability to initialize a large number 
CS) I of qubits into a specific pure state. In some proposed physical implementations, however, this appears to be 

^ ' an extremely difficult task. In 1998, Knill and Laflamme proposed that exponential speedups over classical 

computers could still be possible, even if one can only initialize a single qubit into a pure state, with the 
rest of the qubits in the maximally mixed state [T7]. The complexity class thus defined is called DQCl 
[ (deterministic quantum computation with one clean qubit), or simply "the one clean qubit class." This class 

1^ ' contains several problems for which no efficient classical algorithms are known. The most basic of these is 

I the problem of estimating the trace of a unitary operator. In fact, trace estimation is DQCl-complete; not 

^— I ' only is it in DQCl, but any other problem in DQCl can be reduced to it. 

T-H I Finding natural BQP-complete and DQCl-complete problems is essential to our understanding of the 

^ ■ computational power afforded by quantum computers. Remarkably, BQP-complete problems can be found 

[ in areas of mathematics without a priori close connection to quantum computation. In particular, approxi- 

rS ' mating the Jones polynomial, a famous invariant of links, is a BQP-complete problem [T^ [T?[ \TM fTl [21 [25 ] . 

I The input is an element of the braid group, and the output is an estimate of the Jones polynomial of the 

" ■ ■ so-called plat closure of the braid. Estimating the Jones polynomial of the so-called trace closure of the 

braid is DQCl-complete [25 1 [TB ] . 

Recent work [3J |TS] showed that (the decision version of) approximating certain invariants of 3- manifolds 
is a BQP-complete problem. In this formulation, the input is a so-called Heegaard splitting of a 3-manifold, 
specified as an element of the mapping class group. The output is an estimate of the Turaev-Viro invariant 
of the input manifold. In this article we show that approximating the Turaev-Viro invariant of a 3-manifold 
specified as a mapping torus is a complete problem for the one clean qubit class. In sectionjSl we use the lan- 
guage of Topological Quantum Field Theories (or TQFTs) to explain the mathematical underpinnings of the 
relationship between approximating the Jones polynomial of the plat and trace closures, and approximating 
the Turaev-Viro invariant of Heegaard splittings and mapping tori. 

'Current affiliation: NIST {stephen.jordan&nist. gov) . This work was done at Institute for Quantum Information, Caltech. 
^Institute for Quantum Computation, University of Waterloo (galagic&iqc.ca). 
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Figure 1: A Dehn twist is a 27r rotation about a closed curve. The Dehn twists along the 3g — 1 curves illustrated 
here constitute a standard set of generators for the mapping class group of the genus g surface. 



We assume only a basic understanding of topology and quantum computation. Needed concepts in 
manifold invariants and one clean qubit computation are explained in section [51 Our exposition focuses on 
the Witten-Reshetikhin-Turaev (or WRT) invariant. This is only a matter of convenience, as it is known 
that the Turaev- Viro invariant is equal to the absolute square of the WRT invariant [521 HZl UHl US] • 

2 Background 

2.1 Two- manifolds and three- manifolds 

We begin by setting down a few basic definitions from low-dimensional topology. Recall that an n-manifold is 
a topological spacq^l whose every point has a neighborhood that looks like (i.e., is homeomorphic to) an open 
subset of R" . Simple examples of one-dimensional manifolds include the line M and the circle . Simple 
examples of two-dimensional manifolds include the the plane M^, the sphere S"^, and the torus Ei = x , 
which we can visualize as the surface of a donut. More generally, the surface of a donut with g holes is also 
a two-manifold, which we call the surface of genus g and denote by S^. The genus is a complete invariant 
of surfaces!!: homeomorphic surfaces have the same number of handles (invariance) , and non- homeomorphic 
surfaces have a different number of handles (completeness) . 

The simplest example of a 3-manifold is R'^ itself. A nontrivial example is found by taking the product 
of El with a third circle; the result is the three-dimensional torus = x x . Given a surface S^, 
the cylinder "Eg x [0, 1] is a 3-manifold whose boundary consists of two copies of Eg (specifically, the bottom 
Tig X {0} and the top x {!}.) We can turn the cylinder into a 3-manifold without boundary by choosing 
a homeomorphism / : Eg — > Eg and gluing each point on the top to its image under / on the bottom. The 
result is the mapping torus of /: 

Eg X [0, 1] 

^'^ (x,l)^(/(x),0) • 

For example, choosing g = I and / to be the identity map, we see that Ti^ = T^. A useful example of a 
nontrivial self-homeomorphism of Eg is the so-called Dehn twist. To visualize a Dehn twist, imagine cutting 
the handle of Ei to get a tube, performing a 27r twist on one end of the tube, and then gluing the handle 
back together. In general, a Dehn twist can be performed around any noncontractible closed curve. 

The (homeomorphism class of) the mapping torus Tgj depends only on the isotopy class of /. The 
orientation-preserving self-homeomorphisms of Eg form a group under composition. This group, taken 
modulo isotopy, is called the mapping class group of Eg, and is denoted MCG(g). MCG(g) is generated 
by the Dehn twists about the 3g — 1 canonical curves shown in figure [H Any mapping torus Tgj is thus 
described by a word in the Dehn twist generators of MCG(g). 

^ more precisely, a second-countable Hausdorff space 

^In this work, we implicitly assume that all surfaces are closed, compact, connected and orientable. 
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Figure 2: The dashed hnes indicate a set of cuts that decomposes the surface into two three-punctured-spheres 
("pants"). Dual to this is a trivalent graph called the "spine," in red. The genus two and genus three cases are shown 
here. 




Figure 3: The Fibonacci model's fusion rules allow five labelings of the standard spine of the genus two surface. This 
means that the WRT representation of MCG(2) is five-dimensional. 



2.2 The Witten-Reshetikhin-Turaev invariants 

Recall that the genus is an invariant of surfaces because it assigns the same number to homeomorphic 
surfaces. One can also define invariants of 3-manifolds, although none are as simple and powerful as the 
genus. In the 1990s, Witten, Reshetikhin, and Turaev discovered a family of 3-manifold invariants arising 
from their work in Topological Quantum Field Theory. While these invariants can be defined for arbitrary 
3-manifolds, we only concern ourselves with the special case of mapping tori, where the definitions are 
relatively straightforward. Specifically, the Witten-Reshetikhin-Turaev (WRT) invariant of a mapping torus 
Tgj is equal to the trace of / in a certain projective representation of the mapping class group MCG^g). 
Note that the WRT function is only a topological invariant up to a phase (see [^). In general, the WRT 
invariant is parametrized by a quantum group, such as SU(Af)fe or SO(iV)fc. Although some of our results 
apply more generally, we focus on the case of SO(3)3, sometimes called the Fibonacci model. In this case, 
the description of the representation is particularly simple, and can be understood with no background in 
quantum groups. 

The Fibonacci representation is defined as follows. Any genus g surface (for g > 1) can be cut into 
three-punctured spheres, resulting in a so-called pants decomposition. Dual to such a decomposition is a 
trivalent graph on the surface, called a spine. As illustrated in figure [2l the spine has one vertex for every 
pant in the decomposition. Whenever two pants meet at a puncture, the spine has an edge between the 
corresponding vertices. While a surface admits many spines (and corresponding pants decompositions), we 
call the one shown in figure [2] the standard spine. We label the edges of the standard spine by so-called anyon 
types, with fusion rules enforced at each vertex. For the Fibonacci model, there are only two anyon types: 
and 1, and only one fusion rule: no vertex can have exactly two edges labeled incident on it. The case 
5 = 2 is pictured in figure [H The formal span (over C) of all such labelings associates a finite-dimensional 
vector space to the surface. Different spines yield different bases for this same space. We can move between 
these spines (and the corresponding bases) by means of two "moves," the F-move: 
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and the S-move: 




For the Fibonacci model F^l^^ is as follows 
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with all other values equal to zero or one as dictated by the fusion rules. As one can calculate using the 
prescription described in [3], 5']^. is given in the Fibonacci model by 
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with I? = Y 1 + ^ ^^2^ ) other values of S'jj, equal to zero by the fusion rules. 

The space described above is the underlying vector space for the Fibonacci representation of MCG{g). 
We define this representation in the basis corresponding to the standard spine. Since the mapping class group 
is finitely-generated, it suffices to describe the images of the Dehn twist generators. Any such generator is 
a 2n twist along some canonical curve c from figure [1] It is not hard to check that, by applying at most 
one F-move and one S-move, the standard spine can be adjusted so that c is a cut in the corresponding 
pants decomposition. In this basis, the Dehn twist about c induces a diagonal linear transformation. To 
each labeling of the spine corresponds a basis vector, and this basis vector obtains a phase determined by 
the label on the edge of the spine that intersects c. In the Fibonacci model, edges labeled obtain a phase 
of 1, and edges labeled 1 obtain a phase of e^^^^^. In the standard spine basis, the matrix corresponding to 
the Dehn twist about c is thus simply a product of at most five matrices: at most two of the moves pictured 
above, followed by a diagonal matrix, followed by the inverse moves to return to the original basis. The 
WRT invariant of the mapping torus Tg j is now simply the trace of the Fibonacci representation, evaluated 
at /. 



2.3 One Clean Qubit 

In some proposed implementations of quantum computers, such as nuclear magnetic resonance (NMR) the 
most difficult task is initializing qubits into a pure state. In 1998, Knill and Lafiamme proposed that 
exponential speedups over classical computation might be possible without pure state initialization. To 
mathematically investigate this possibility, they introduced the one clean qubit model |17j . In this model. 
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Figure 4: By repeating this one clean qubit computation, and recording the fraction of outcomes, one estimates 
the real part of Tr[(7]/2". Similarly, by initializing the clean qubit to -^(10) — i |1)), one obtains po — ^ + ^"'(^I'-^D , 



one is given an initial state p with n qubits in the maximally mixed state, and one qubit in the pure state 
|0). 

P = |o> (0| ® ^ 

One then applies any quantum circuit of poly(n) gates to this state, and measures the first qubit in the com- 
putational basis. Computational problems are solved by performing polynomially many such experiments, 
each starting with the initial state p, and recording the output statistics. The class of decision problems 
solvable with bounded probability of error using this procedure is called DQCl. 

DQCl contains several computational problems not known to be solvable in polynomial time on classical 
computers. Most fundamentally, given a description of a quantum circuit of T gates implementing the unitary 
transformation U onn qubits, a one clean qubit computer can estimate the normalized trace to within ±e 
in time 0(T/e^) by means of the circuit shown in figured Furthermore, this problem of estimating the trace 
of a quantum circuit is DQCl-hard [T71 [231 [IS]. Efficient one clean qubit algorithms have been discovered 
for estimating certain quadratically signed weight enumerators [T^ and estimating certain Jones [3S] and 
HOMFLY [12] polynomials. A version of the Jones polynomial problem is DQCl-complete [3S], and has been 
demonstrated experimentally with NMR |221 120] . A certain problem of approximating partition functions 
for quantum systems is also DQCl-hard [6]. 

In many ways, it is surprising that one clean qubit computers can do any nontrivial computations at 
all. If all n -f 1 qubits were maximally mixed, the resulting state would be invariant under all unitaries. 
Furthermore, DQCl computations involve very little entanglement [7] [HI [SJ [TUl [HJ Ambainis et al. 

give an impossibility proof against a certain natural approach to simulating standard quantum computers 
using one clean qubit computers, and on the other hand show that one clean qubit computers can efficiently 
simulate classical logarithmic depth (NCI) computations [4]. 

The DQCl complexity class is robust against a variety of modifications to the computational model. The 
class of computational problems solvable in polynomial time with up to logarithmically many clean qubits 
is the same as that solvable in polynomial time with one clean qubit |25) . If the clean qubit is not pure, 
but has l/poly(n) polarization, the set of efficiently solvable problems also remains DQCl [17]. As shown 
in appendix [3 the one clean qudit model on d-dimensional qudits is equivalent in power to the one clean 
qubit model, for any constant d. 

3 Algorithm 

In this section we construct an efficient one clean qubit algorithm for approximating the Fibonacci WRT 
invariant of a mapping torus. Generalizing to other tensor categories such as SU(A'')fe and SO{N)k is 
straightforward. The main idea of the algorithm is, given a word w in the Dehn twist generators of MCG^g), 
to find a quantum circuit of poly(w,g) gates on poly{g) qubits whose trace is equal to the WRT invariant 
of the 3-manifold Tg^^. This trace can then be approximated by means of the circuit in figure [4] For this 
purpose, we encode the allowed labelings of a spine of Eg into qubits, and then construct a quantum circuit 
implementing the Fibonacci representation of MCG((7) on this encoding. The most obvious encoding would 
be to directly assign one qubit to store the particle type for each edge of the spine. However, a one clean qubit 
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Figure 5: We number the edges of the standard spine from left to right, with ambiguities resolved by ordering from 
top to bottom. 



computer yields the normalized trace over all 2" bitstrings, of which only an exponentially small fraction 
represent valid spine labelings in this encoding. 

We instead construct a many-to-one map (/9 : {0, — {valid labelings} with /? = 0(log |g|) such 
that the preimage of each spine-labeling consists of approximately the same number of bitstrings. That is, 
is approximately independent of x. Thus, the normalized trace of the Fibonacci representation of 
w S MCG(g) acting on the (^-encoded labelings of the spine of Eg is approximately equal to WRT(Tg,TO). 
We construct if following a method introduced in jJL6j. We assign a register of /3 = 0(log \g\) qubits to each 
edge of the spine. The bitstring contained in register i is interpreted as an integer Q < Xi < 2^ — 1. We then 
assign a threshold Ti so that Xi < Ti encodes a zero label on edge z, and Xi > Ti encodes a one label. By 
carefully choosing the thresholds Ti, . . . ,T3g_3 we ensure that is approximately independent of x. 

Number the edges of the spine from one to 2>g — 3, left to right and top to bottom, as illustrated in figure 
[3 Let Si, . . . , S3g-3 G {0, l}3s-3 be the labels of these edges. The uniform probability distribution over all 
fusion-consistent labelings of the spine induces a probability distribution Pg{si . . . ,533-3) over {0, \}^^~^, 
with zero probability for strings that violate fusion rules, and uniform probability for the rest. For the 
genus-g standard spine, we define pg{si\si, . . . , st-i) to be the conditional probability that label i takes the 
value Si given that labels 1 through i — 1 take the values si, . . . , Si_i. For a register representing a label Si 
we choose the threshold dependent on the values of si, . . . , si-i according to 

T,{g- si, . . . , = [2^p<,(0|si, . . . , s,_i)J . (1) 

One can see that this choice ensures that a uniformly selected assignment of bitstrings to the registers yields a 
uniform distribution over fusion-consistent labelings, up to the errors induced by rounding. Hence, |(/9^^(a;)| 
is approximately independent of x. More precisely, let 

Pg(0|si,...,Si_i) = Ti{g]Si,...,Si-i)/2^ 
Pg(l|si,...,Sj_i) = 1 -pg(0|si,...,Si_i) 

2''^^''"^^Pg(s3g-3|si, ■ • • , S3ff-4) X _Pg(s3g-4|si, . . . , SSg-s) X ... X p{si) 
2^(39-3) (p^(s3^_3|s,, . . . , S3^_4) ± 0(2-^)) X ... X (pg(si) ± 0(2-^)) 

Pg(si,...,S3ff-3)±0(52-'5). 

Oilogg). Furthermore, by the locality of the fusion rules, Pg{si\si, . . . , Si-i) 
, Si„3. We may thus write 

Pg(Si|si, . . . , Si_i) = Pg(si,Si_i,Si_2;«) 

Ti(g; Si, . . . , Si_i) = T^{g]Si,s.i-i,Si-2)- (2) 

As illustrated in figure [6l the Fibonacci representation of a Dehn twist from the standard generating set 
is a unitary transformation acting on at most five spine labels. Because the encoding ip is many-to-one, the 
unitary transformation on these spine labels does not uniquely define a unitary operation on the bitstrings 
encoding them. We say that a pair of spine-labelings is connected if the Fibonacci representation of a Dehn 
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6 




acts on a 



a 



e 




acts on a, b, c, d, e 



b 



d 



Figure 6: The Fibonacci representation of a Dehn twist (shown as a dashed line) from the standard generating set is 
a unitary transformation acting on at most five spine labels. 



twist from the standard set of generators has a nonzero matrix element between them. By choosing a bijection 
bx,y between the encodings of each pair of connected spin-labelings we define a unitary transformation on 
the encodings: if the matrix element between labeling x and y is px.y then, 



is a corresponding unitary representation on the encodings. Our choice of bijections does not matter. We 
may for concreteness match bitstrings by lexicographic ordering. One can verify that U is a, direct sum of 
many copies of the Fibonacci representation p. (The rounding involved in ([T|) introduces a minor technical 
complication, whose resolution may be found in [16].) 

For any of the standard Dehn twist generators, Uij acts on at most 5/? qubits, which encode the spine- 
labels on which p acts. The matrix elements by which U acts on these qubits depends on the corresponding 
thresholds. By these depend on at most two additional registers of qubits, which encode the two spine 
labels to the left of those being acted upon. Thus, for any of the standard Dehn twist generators, Uij is a 
controlled unitary acting on at most 5/3 target qubits and 2/3 control qubits. Recalling that (3 — 0(log j^l), 
we can apply the standard construction from section 4.5 of [21^ to implement this unitary transformation 
with polyd^l) quantum gates, provided each matrix element of Ui,j can be computed efficiently. By (jSj, one 
sees that the only potentially difficult part of computing the matrix elements of [3] is the computation of the 
thresholds. An efficient classical algorithm for this task is given in appendix [Cl 

4 Hardness 

In this section we prove that the problem of estimating the normalized WRT Fibonacci invariant of a mapping 
torus, given by a polynomial-length word in the standard Dehn twist generators of the mapping class group, 
to within ±e is DQCl-hard for e < 1/3900. Generalizing our hardness proof beyond the Fibonacci model 
seems less straightforward than generalizing our algorithm. However, we consider it likely to be possible. 
Extending hardness to larger values of e we leave as an open problem. To prove hardness, we reduce from 
the problem of estimating the absolute value of the normalized trace of a quantum circuit. A proof of the 
hardness of absolute trace estimation is given in appendix IbI We thus require an efhcient procedure that, 
given a description of a quantum circuit for implementing a unitary U , outputs a description of a mapping 
torus (i.e., a word in the Dehn twist generators) whose WRT invariant is close to the trace of U. It turns 
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out to be convenient to suppose that /7 is a quantum circuit acting on a collection of 5-dimensional qudits 
("qupents"). As shown in appendix |X1 this makes no difference: the one-clean-qubit model is equivalent to 
the one-clean-qupent model. 

Let U he a quantum circuit of G gates acting on n qupents arranged in a line. Without loss of generality, 
we may assume that each gate acts either on a single qupent or a pair of neighboring qupents. To prove 
hardness, we first define a many-to-one encoding ■0 : S^n — ^ {Oj 1: 2, 3, 4}", where S^n is the set of fusion- 
consistent labelings of the standard spine of the surface of genus 3n. We divide the genus-Sn surface into n 
segments, each having three handles. The number of fusion-consistent labelings for a genus-three segment 
with two punctures depends on the labels on the incoming and outgoing edges, as shown below. 



-0"-0""0~^ labelings 
^ ^0"CH3^ ^ '^^ labehngs 
^^0""CH3~° has 20 labelings 

1 -^^^^H^^^^H^^^ 1 has 35 labelings 



In all cases, the number of fusion-consistent labelings is a multiple of five. Thus, in every case a qupent 
can be encoded in the space of labelings, together with a "gauge" qudit, whose value we ignore, which has 
dimension 3, 4, or 7, depending on the labels of the incoming edges. Thus the size of the preimage 

of any z G {0,1,2,3,4}" is exactly independent of z. Given any unitary U acting on n qupents, there 
corresponds a unitary acting on the span of S'3„ which acts as U on the encoded qupents, and as the identity 
on the gauge qudits. We call this the ?/)-encoding of U. 

As shown in [M] , the Fibonacci representation of the mapping class group of the genus g > 1 surface is 
dense in the corresponding unitary group, modulo phase. Thus, given any unitary operation on n qupents, 
we can find a sequence of Dehn twists which approximates its '(/'-encoding arbitrarily closely. The trace of 
the i/i-encoding is thus equal to the trace of the original quantum circuit, up to a phase. The remaining 
question is whether this reduction can be done efRciently. 

Cutting the genus-3n surface into n equal segments yields n — 2 genus-3 doubly-punctured surfaces, and 
two genus-3 singly-punctured surfaces, as shown below. 

x(n- 2) 



^<:>' vcy 



VO' vcv VZV vc>' ^z^- 



One can pants-decompose a punctured surface, thereby associating the surface to a spine. The spine has 
one "external" edge for each puncture, which attaches to the rest of the spine at only one vertex. Upon 
labeling the spine, we can associate the label of any external edge with the corresponding puncture. The 
Fibonacci representation may then be extended to the label-preserving mapping class group of the punctured 
surface. This group includes all the standard Dehn twists, together with braiding of punctures with the other 
punctures of the same label. In the Fibonacci representation, braiding of zero-labeled punctures has no effect, 
thus a zero-labeled puncture is equivalent to the absence of a puncture. 

Theorem 6.2 of [l4, states that for any fixed labels on the punctures, the Fibonacci representation of the 
label-preserving mapping class group of the r-punctured genus-y surface is dense in the corresponding unitary 
group modulo phase, provided g + r > 1. Thus, given any one-qupent gate, the Solovay-Kitaev theorem [21] 
efficiently yields a sequence of Dehn twists and braid moves on the corresponding genus-3 singly-punctured 
or doubly-punctured surface, whose Fibonacci representation approximates the -(/-encoded gate arbitrarily 
closely. Similarly, one efficiently approximates two-qupent gates by moves on genus-6 surfaces with one or 
two punctures. 

We must modify the above construction so as not to use any braiding of punctures. On the leftmost or 
rightmost qupents there is no problem; the corresponding surfaces have only one puncture, and therefore 
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Figure 7: The problems of estimating the Jones polynomial of the plat closure of a braid and the Turaev-Viro 
invariant of a Heegaard splitting (left) are BQP-complete. The problems of estimating the Jones polynomial of the 
trace closure of a braid and the Turaev-Viro invariant of a mapping torus (right) are DQCl-complete. These situations 
are fundamentally analogous, as discussed in section [5] We stress that the manifold figures are illustrations of the 
topological ideas behind this analogy, and are not correct two-dimensional projections of the manifolds themselves. 
In particular, after gluing, the two manifolds shown do not in reality have any boundaries. 



theorem 6.2 implies density without using any braiding moves. Similarly, on any of the central surfaces, 
theorem 6.2 implies density without using any braiding moves if at least one of the punctures has a zero 
label. We can ensure this prior to the application of any given gate by adapting the "inchworm" technique 
from as described in appendix iDl In this method, we bring a pair of zero labels adjacent to the target 
segment, then implement the desired gate there, and carry the zeros to the segment where the next gate 
is to be implemented. At the end, we return these zeroes to their original location among the leftmost six 
handles. As discussed in appendix[Dl the inchworm construction entails some overhead in e, which gives rise 
to the value 1/3900. 

In the above construction, we need density on two-punctured segments in which one puncture is guar- 
anteed to be labeled zero, and the other puncture has unknown label. Theorem 6.2 of |14| implies density 
separately in the subspace in which the other label is zero and in which the other label is one. Because these 
subspaces have different dimension (20 and 15, respectively) we may apply the decoupling lemma from [T], 
which shows that a sequence of Dehn twists can be found to approximate arbitrary pairs of independent 
unitaries on these two subspaces, as desired. 



5 Analogy with Jones polynomials 

In this paper we have shown that estimating the Turaev-Viro invariant of a mapping torus in the Fibonacci 
model is DQCl-complete. In [3], it was shown that estimating the Turaev-Viro invariant of a general 3- 
manifold presented as a Heegaard splitting is BQP-complete. Similarly, estimating the Jones polynomial of 
the trace closure of a braid is DQCl-complete [25l [16], while estimating the Jones polynomial of the plat 
closure of a braid is BQP-complete [H [TJ [521 [T31 [H] . This suggests a relationship between trace closures and 
mapping tori on one hand, and between plat closures and Heegaard splittings on the other. Indeed, such a 
relationship can be understood in the framework of axiomatic topological quantum field theory, and suggests 
further generalizations to, for instance, topological invariants of higher dimensional manifolds. 

A topological quantum field theory can be axiomatized as a functor T from the category of cobordisms 
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Figure 8: M can be viewed as a two-manifold with two boundaries: a circle and a pair of circles. The TQFT 
associates a Hilbert space T{C2) to the pair of circles, a Hilbert space T(Ci) to the circle, and a linear transformation 
T(M) : T(C2) r(Ci) to M. Similarly, T(N) is a linear transformation from T(Ci) to itself. If we glue together M 
and A'^ along the circle as shown, we obtain a manifold MN with boundaries C2 and Ci. The corresponding linear 
transformation is T{NM) = T{N) o T{M). 



between n-manifolds to the category of linear transformations between vector spaces [51 [35]. That is, to each 
n-manifold the TQFT associates a vector space, and to any (n + l)-manifold whose boundary is the union 
of two disjoint n-manifolds the TQFT associates a linear transformation between the two associated vector 
spaces. The functorial property means that gluing together two cobordisms and then applying T yields the 
same Hnear transformation that is obtained by applying T to each of the two cobordisms and then composing 
the resulting linear transformations; see figure [Sj A TQFT maps the empty n-manifold to the base field, 
which for the examples we consider is always C. Hence, for M a manifold whose boundary dM has a single 
connected component, T{M) is a map either from C to the vector space T{dAI), that is, a vector in T{dM), 
or a map from T{dM) to C, that is, a dual vector. The choice between these two possibilities is determined 
by the orientation of the cobordism. 

Recall that the genus-g handlebody is the 3-nianifold whose boundary is the genus-g surface Eg. For 
example, the genus-1 handlebody is simply the solid donut. After assigning an orientation, we may think 
of a handlebody as a cobordism from the empty manifold to Sg, or as a cobordism from Eg to the empty 
manifold. Hence, in the TQFT framework, genus-g handlebodies are associated to vectors or dual vectors. 
We denote these as j^'g) and (-(Agl, respectively. These vectors live in the Hilbert space which the TQFT 
associates to Eg. In the case of the Fibonacci model, this is precisely the vector space defined in Section [2?2] 

In the Fibonacci model, a cobordism from a surface to itself is mapped to a unitary linear transformation 
U on the associated Hilbert spac^. If the surface is Eg, then we may "cap" the cobordism with handlebodies 
on both ends. The resulting 3-manifold has no boundary, and thus corresponds to a linear map from C to 
itself, i.e., a complex number. In this case, this number is the matrix element {'ipg\ U \'lpg)^ as illustrated in 
figure [7] The problem of estimating a matrix element of the unitary transformation induced by a quantum 
circuit is BQP-complete, and this fact underlies the BQP-completeness proof for the Turaev-Viro invariant 
of Heegaard splittings in [3]. Instead of "capping" the two ends of the cobordism with handlebodies, we 
could have simply glued the two ends together, resulting in a mapping torus. This is again a 3-manifold 
without boundary, which thus also corresponds to a complex number. In a TQFT, gluing the two ends of a 
cobordism corresponds to contracting the two indices of the linear transformation. In other words, instead 
of a single matrix entry, we now obtain the trace of U . Finding the trace of the unitary transformation 
induced by a quantum circuit is DQCl-complete, and this fact underlies the DQCl-completeness proof for 
the Turaev-Viro invariant of mapping tori given in this paper. 

The situation regarding Jones polynomials is directly analogous. A TQFT gives us a unitary representa- 
tion of the braid group. Gluing the two ends of a braid together (i.e., taking the trace closure), as illustrated 
on the righthand side of figure[7l corresponds to taking the trace of the unitary and yields a DQCl-complete 
problem. Caps correspond to vectors and dual vectors depending on orientation, hence capping a braid 

^We may think of the cobordism as describing a sort of spacetime evolution, while the unitary transformation describes the 
corresponding quantum time evolution. Indeed, this was one of the central motivating ideas behind the development of TQFTs. 
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(taking the plat closure, as illustrated on the lefthand side of figure [7]) yields a matrix element of the asso- 
ciated unitary transformation, and corresponds to a BQP-complete problem. The analogy can be tightened 
further by noting that the braid group i?„ is simply the mapping class group of the surface of genus zero 
and n + 1 punctures (that is, the n-punctured disk) , whereas in the case of 3- manifold invariants we consider 
the mapping class group of the genus-g surface with no punctures. On the other hand, it is worth bearing 
in mind that the notion of equivalence captured by the Jones polynomial is ambient isotopy, in contrast to 
the Turaev-Viro and WRT invariants, which capture homeomorphism. 

The analogy presented here naturally suggests an extension of BQP-completeness and DQCl-completeness 
results to rt- manifold invariants arising from TQFTs at higher n. More generally, one could attempt to isolate 
a property of pairs, consisting of a group G and one of its representations U, such that estimating matrix 
entries of U is BQP-complete while estimating the trace of U is DQCl-complete. Perhaps one could find a 
general theorem encompassing many such results. We leave this as an open problem. 
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A Equivalence Between One Clean Qudit Models 

Given a quantum circuit on a-dimensional qudits we wish to construct a quantum circuit on 5-dimensional 
qudits that has the same trace. If 6 = ca for some integer c then this is easy. We just consider each b- 
dimensional qudit to be an a-dimensional qudit plus a c-dimensional "gauge" qudit that we ignore. Similarly, 
if b'^ = ca for some integers d, c then we can treat c?-tuples of 5-dimensional qubits as an a-dimensional qudit 
plus a c-dimensional gauge qudit. For these encodings, the encoded circuit is easy to construct gate by gate. 
Given a gate acting on n a-dimensional qudits, we can write down a unitary acting on dn 5-dimensional qudits 
equal to the original gate tensored with the c-dimensional identity on the gauge system. This dn-dimensional 
gate can be exactly decomposed into a product of 0(6^*') 2-qudit gates using the standard construction 
from section 4.5 of 21 . Because d and n are constants, this is sufficiently efficient. The normalized trace of 
the encoded circuit is exactly equal to the normalized trace of the original circuit. 

The harder case is when there do not exist integers c and d such that b'^ = ca. In this case we find c, d € Z 
such that b"^ ~ ca. Specifically, suppose we achieve 

for some <C 1. Then we can encode one a-dimensional qudit plus a c-dimensional gauge qudit into d 
&-dimensional qudits with a few (namely 56'') noncoding states left over. We can define our encoded gates 
to act as the identity on these noncoding states. If we make sure the noncoding states are a small fraction 
of all 6''" states, the normalized trace of the encoded circuit will approximately match the normalized trace 
of the original circuit. 

Let Ua be the original unitary acting on n a-dimensional qudits and let Ub be the unitary acting on dn 
^-dimensional qudits, in which we encode Ua as described above. Then, Ub acts on 6''" states, of which (ca)" 
encode states of the original circuit, 

Trpb] _ c"Tr[;7a] + (fo''" - (ca)") 

^dn ^dn 
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The magnitude of the discrepancy A between the normahzed traces of Ub and Ua is thus 



< 



< 



c"Tr[C/,] + (fo"^" - (ca)") Tr[C/J 



- 1 



1 



- 1 



a" 
Tr[C/J 



1 - 



2|(l-5)"-l| 



1 - 



Thus if 



we have, for small e, 



e 

n 



lim A 



2 e" 



1 



2e. 



(5) 



(6) 



Comparing (jl)), ([5|), ([6]), we see that in the limit of large n and small e, in order to achieve error upper 
bounded by A it suffices to obtain 

b'^-ca A 
I? - 2n' 



For given b, d, a there always exists an integer c such that - 
large that 

a ^ A 
¥ ~ 2n' 



c < a. So we just need to choose d sufficiently 



Equivalently, 



d > log,, 



2na 
"A~ 



dk 



A fc-qudit gate from Ua thus gets encoded as a d/c-qudit gate in Ub- This encoded gate acts on a b' 
dimensional space. We have just shown that it suffices to choose d — log^ (^s^)- Thus the encoded /c-qudit 

gate acts on a (222^)'°-dimensional space. Using the construction from section 4.5 of |3T], we can implement 
an arbitrary /^-dimensional unitary exactly with 0{D^) 2-qudit gates. Thus each fc-qudit gate in Ua gets 



encoded by O 



elementary gates in Ub- By gate universality, we can assume fc < 2, so our encoding 

has an overhead quartic in n and 1/A. This is perhaps not very efficient, but is nevertheless polynomial, 
and thus sufhces to prove the equivalence of DQCl defined with qudits of any constant dimension. 



B Estimating the Absolute Trace is DQCl-hard 

In this section we slightly adapt the proof from to show that estimating the absolute value of the trace of 
a quantum circuit to within ±1/24 is a DQCl-complete problem. Consider an arbitrary DQCl computation. 
We start with the state |0) (0| iX) apply an arbitrary quantum circuit J7, and then measure the first qubit 
in the |0) , |1) basis. Changing the initial state of the pure qubit, or changing the measurement basis does 
not add generality, as these changes can be subsumed into U- The probability of measurement outcome |0) 
is 

po = Tr [(|0) (0| ® t)Um (0| ® l/2")C/t] . (7) 
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Let U' be the unitary implemented by the fohowing quantum circuit on n + 2 qubits. 
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Thus, po = 2^^^, as one can see by writing out the trace as a sum over diagonal matrix elements in the 

computational basis. Because po is real it is also true that po — ^-^pr- Hence estimating the absolute value 
of the normalized trace of quantum circuits to suffices to predict the outcome of any DQCl experiment. 

As is standard in the complexity theory of probabilistic computation, "yes" instances of DQCl are defined 
to have acceptance probability 2/3 and "no" instances are defined to have acceptance probability 1/3. Thus, 
deciding DQCl is equivalent to estimating the normalized trace of a quantum circuit to within ±1/6. The 
reduction here has a factor of four overhead in normalization, thus estimating the absolute trace to within 
±1/24 is DQCl-complete. 



C Efficiently Computing Thresholds 

Consider the standard spine of the genus-g surface, numbered as in figure [5] Suppose edges 1 through i 
have been labeled in a fusion-consistent manner with anyon types si, . . . , s^. We wish to compute how many 
completions there are to this partial labelling. That is, we wish to compute the number of fusion-consistent 
strings of 3g — 3 labels, whose first i labels are given by si, . . . , s^. 

Denote the horizontal edges of the standard spine from right to left by ei, 62, . . ., as shown below. 



62 



Let Z^^'^ be the number of completions in which the rightmost labeled edge is Cfe and has label h e {0, 1}. 
One sees that Zq^'' = 2, and zj^'' = 1, by the following enumeration of fusion-consistent diagrams. 



Furthermore, we have the recurrence relations 4"+^^ = 2^,$"' + and 
the following enumeration of fusion-consistent diagrams. 



= 3Z 



+ ^r'\by 
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Solving these recurrence relations yields 



^0 



2 1 
1 3 



n-l 



The other two cases — completions starting with an upper curved edge, or a lower curved edge — can be 
solved similarly. The n"^ power of a matrix may be computed using O(logn) operations, thus calculating 
the number of completions for any i in 0{\.ogg) steps. The corresponding thresholds are then immediately 
obtained by taking ratios of these. 



D Inchworm 

Suppose the spine-labeling contains a segment of the following form. 




(8) 

Here a can be 1 or 0. We call this configuration the inchworm. We may regard the right instance of 
1 1 

— as its head, and the left instance as its tail. We next show a sequence of two reversible 
1 1 

operations by which we can move the inchworm one handle rightward. In the first step the head moves one 
handle to the right, leaving the tail in place, and in the second step, the tail catches up, hence the name 
"inchworm." 




Examination of the above diagram shows that if the fusion rules arc obeyed in the initial configiiration, they 
are also obeyed in the intermediate and final configurations. Furthermore, both steps are reversible (i.e. 
information preserving). Thus, they may be written as permutation matrices acting on the space of allowed 
configurations, and are therefore unitary. The first unitary transformation can be implemented by local 
Dehn twists, because the zero in the tail of the inchworm implies density of the Fibonacci representation 
on the segment to the right of it. The second unitary transformation can be implemented by local Dehn 
twists because the zero in the head of the inchworm implies density on the segment to the left of it. (In both 
steps, we are applying density to the twice-punctured genus-4 surface with one puncture labeled zero. There 
are 75 labelings in which the other puncture is labeled one and 50 labelings in which the other puncture is 
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labeled zero. Thus, the decoupling lemma of [T] implies density jointly on these two subspaces.) Repeating 
this process and its reverse, we may bring the inchworm to any location within the spine. 

To use the inchworm construction, we need to ensure that a segment of the form ([5]) exists in the 
first place. We may do this by implementing a reversible operation on the leftmost six handles, so that 
if the configuration ([5]) is absent, the matrix is strictly off-diagonal, and does not contribute to the trace. 
Specifically, we consider the leftmost two handles to be an ancilla system, and the next four handles to be 
the starting location of the inchworm. If these four handles do not take the form ([8]) we cyclically permute 
the (five) basis states of the ancilla system. Because this is done on the leftmost six handles, the segment is 
only singly-punctured, and thus theorem 6.2 of |14| implies density without braiding. 

The noncontributing labelings decrease the normalized WRT by a constant factor, which correspondingly 
necessitates decreasing the precision parameter e by the same factor. More precisely, in the Fibonacci model, 
there are 325 fusion-consistent labelings for the spine of the genus-four doubly-punctured surface. Among 
these, there are two inchworm configurations {a — and a = 1). Compounding this 2/325 normalization 
cost with the precision e = 1/24 obtained in appendix |B] for DQCl-hardness of absolute trace, we find that 
estimating the normalized WRT invariant to within ±1/3900 is DQCl-hard. 

As an aside, we note that the inchworm construction here is simpler than that in [33], in the following 
sense. The inchworm construction of [2 5) involved reversible operations on logarithmically large regions. 
Although the density theorems imply that arbitrary reversible operations can be implemented on these 
regions, they do not imply that the decomposition into local moves is efficient. Rather this had to be 
explicitly proven in appendix D of |25| . In contrast the inchworm construction here involves reversible 
operations only on 0(1) handles, thus no question of efficiency arises. 
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